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ABSTRACT 

We present composite UV spectra for low redshift Type 1 AGN binned to ex- 
ploit the information content of the Eigenvector 1 (El) parameter space. Com- 
posite spectra show high enough S/N and spectral resolution to permit a de- 
composition of the CivA1549 line profile - one of the strongest high-ionization 
lines (HIL) and fundamental in defining El space. The simplest CivA1549 de- 
composition into narrow (NLR), broad (BLR) and very broad (VBLR) compo- 
nents suggests that different components have an analog in H/3 with two major 
exceptions. VBLR emission is seen only in population B (FWHM(H/3bc)>4000 
km s _1 ) sources. A blue shifted/asymmetric BLR component is seen only in pop- 
ulation A (FWHM(H ( 5 B c)<4000km s' 1 ) HIL such as CivA1549. The blueshifted 
component is thought to arise in a high ionization wind or outflow. Our analysis 
suggests that such a wind can only be produced in population A (almost all radio- 
quiet) sources where the accretion rate is relatively high. We propose a model to 
account for several differences between low- and high-ionization line profiles. Part 
of the broad line emission is attributed to a self-gravitating/fragmented region 
in an accretion disk. An inner optically thick geometrically thin region of the 
flow may give rise to a wind/outflow and produce the blueshifted HIL spectrum 



department of Physics and Astronomy, University of Alabama, Tuscaloosa, AL 35487, USA; gia- 
como@merlot.astr.ua.edu 

institute of Astronomy, Sofia 1784, Bulgaria; bachevr@astro.bas.bg 

3 INAF, Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, 1-35122 Padova, Italy; 
marziani@pd.astro.it; calvani@pd.astro.it 

4 National Astronomical Observatory Rozhcn, POB 136, Smoljan, BG-4700, Bulgaria 

5 Instituto de Astronomfa, UNAM, Apdo. Postal 70-264, 04510 Mexico D.F., Mexico; debo- 
rah@ ast roscu . unam . mx 



- 2 - 



in Pop. A sources. The fragmented region may produce all or most of the broad 
line emission in population B, which contains radio-quiet and the majority of 
radio-loud quasars. Comparison between broad UV lines in radio-loud (RL) and 
radio-quiet (RQ) sources in a single well populated El parameter space bin (Bl) 
shows few significant differences. Clear evidence is found for a significant NLR 
CIV component in most radio-loud sources. The BLR/VBLR similarity in bin 
Bl provides circumstantial evidence in favor of black hole (BH) spin, rather than 
BH mass or accretion rate, as a key trigger in determining whether an object will 
be RL or RQ. We find a ten-fold decrease in EW CivA1549 with Eddington ratio 
(decreasing from «1 to ~0.01) while NvA1240 shows no change. These trends 
suggest a luminosity-independent "Baldwin effect" where the physical driver may 
be the Eddington ratio. 

Subject headings: quasars: emission lines - quasars: general 

1. Introduction 

Understanding the geometry and kinematics of the central regions in AGN is one of the 
major goals of twenty-first century astronomy. Yet little is known beyond the necessity of 
invoking accretion onto a supermassive object in order to account for the enormous source 
luminosity. Given our inability to resolve the central regions of even the nearest AGN we are 
left with three approaches towards achieving our goal: (1) ideology: sophisticated models 
driven entirely by theories describing black hole and associated accretion disk properties 
(e.g., Janiuk, Zycki, & Czerny 2000; Gammie, Shapiro, & McKinney 2004), (2) analogy: 
comparison with better-resolved galactic objects (e.g. interacting binary stars; Zamanov & 
Marziani 2002; Mirabel 2004; Maccarone, Gallo, & Fender 2003) and (3) empiricism: detailed 
studies of phenomenology in large and diverse samples of AGN (see Sulentic, Marziani & 
Dultzin-Hacyan 2000a for references until late 1999; Marziani et al. 2003a, hereafter M03; 
Dietrich et al. 2002; Kuraszkiewicz et al. 2004 for examples of more recent works). The 
last approach fell on hard times until the end of the twentieth century in part due to: 
(i) the apparent intractability of AGN phenomenology, (ii) the failure of the best nebular 
physics models to reproduce the observed spectra, and (iii) the low quality of so much of 
the spectroscopic data. Thus we now know more about the morphology of the extended 
radio and optical emission around quasars than the underlying geometry/physics that drives 
them. One positive development involves the application of reverberation mapping towards 
estimation of the broad line region size (see Peterson 1993 for a thorough review; Home et 
al., 2004 for up-to-date expectations). The relevant questions are how to interpret these data 
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and, more importantly, how to select a representative sample of sources for a reverberation 
mapping study. Another positive development involves growing databases of AGN spectra 
(most notably from the Sloan Digital Sky Survey (SDSS); Richards et al. 2002; Schneider 
et al. 2002). Although great part of the individual data still suffers from S/N limitations, 
average or composite spectra of many objects can significantly improve the situation (Vanden 
Berk et al. 2001; Sulentic et al. 2002). One of the most important questions is how 
to bin the data. Clearly unrestricted composites have little value because they hide the 
intrinsic dispersion in key parameters which provides clues to the underlying physics (e.g., 
Francis et al. 1991). In a previous paper (Sulentic et al. 2002: hereafter Paper I) we 
argue that composite spectra binned within the Eigenvector 1 parameter space provide the 
most physically useful results. Paper I presented composites of the H/3 spectral region and 
this paper presents complementary composites for the UV spectral region and especially for 
CivA1549. 

We have identified a four dimensional parameter space that overcomes some of the prob- 
lems mentioned above and that offers a new perspective on AGN phenomenology (Sulentic 
et al. 2000a,b). The Eigenvector 1 (El) parameter space provides optimal discrimina- 
tion/unification of broad-line AGN classes (Boroson & Green 1992; Marziani et al. 2001, 
hereafter M01; Boroson 2002). The adopted name reflects recognition of the first two pa- 
rameters in the pioneering work of Boroson & Green (1992). El as we now define it involves 
measures of: 

1. full width at half maximum of low ionization (LIL) broad emission lines [FWHM(H/3 B c)] , 

2. equivalent width ratio of optical Fell and H/3 BC : R F eii=EW(Fell A4570 blend) /EW(H/3 B c), 

3. soft X-ray photon index (r so f t ), 

4. profile centroid velocity displacement (at FWHM) for CivA1549. 

In simplest terms these parameters can be said to measure respectively: (1) velocity 
dispersion of LIL emitting gas, (2) relative strength of LIL features thought to arise in the 
same region, (3) strength of a soft X-ray component and (4) amplitude of systematic motions 
in the high-ionization line (HIL) emitting gas. In less simple (i.e. more model dependent) 
terms we likely have: (i) two variables mainly sensitive to accretion rate (FWHM(H/3bc) 
and r s : Nicastro 2000; M01; Pounds, Done & Osborne 1995; Czerny et al. 2004), (ii) three 
variables sensitive to source orientation (all with the possible exclusion of Ts; e.g., Jackson 
k Browne 1991; Wills et al. 1995; Marziani et al. 1996; M01; Rokaki et al. 2003; Sulentic et 
al. 2003). RfcII is sensitive to nebular physics (electron density n c and column density N c ). 
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An important addition that reflects the striking changes in emission line properties 
across El involves introduction of the population A and B concept. At present we arbitrarily 
define Pop. A and B quasars as those with FWHM H/3bc l ess than or greater than 4000 
km s _1 . It is unclear whether the A-B concept reflects two distinct quasar populations 
or simply the difference in mean properties of sources at opposite ends of the El "main 
sequence". The "elbow" shape of the sequence, and optical/X-ray indications (Sulentic et al. 
2000a,b; Komossa & Meerschweinchen 2000) for a "zone of avoidance" near FWHM(H/3 B c) ~ 
4000 km s" 1 , provide some evidence in support of the former interpretation. If the distinction 
is real then "A" is a largely RQ population (including narrow line Seyfert 1 sources as 
an extremum), while "B" contains a mix of broader line RQ plus most RL sources. The 
population A-B distinction will prove useful in this study and CivA1549 will add support 
for the concept (§3.2 and §4.3). 

If the El concept lives up to initial promise, it may be the strongest and most encom- 
passing type-1 AGN unification yet discovered. This possibility warrants strong tests of its 
strength and generality. So far most of our effort has centered on the "optical plane" of El 
(M01; Zamanov et al. 2002; Boroson 2002; Marziani et al. 2003a,b; Sulentic et al. 2003). 
This focus has been dictated by the available X-ray and UV archival data that provide the 
other El parameters. In the next section we present composite spectra for CivA1549 (and 
other UV lines) in the El context which complements work already done for E_(3. 



2. Generating Composite UV Spectra 

A total of about 700 reasonable quality (S/N 1 > 4-5) HST spectra for n = 141 different 
AGN (quasars or bright Seyfert l's; n=61 are radio- loud with Rk > 100, where Rk is the 
specific flux ratio between 4400 A and 6 cm) have been extracted from the Hubble archive. 
The majority of the spectra (83 %) were obtained with the FOS camera and the rest with 
STIS. The FOS/BL data (40 %) were recalibrated with the latest IRAF package - STPOA, 
which takes into account some systematic errors in the wavelength calibration. No such 
recalibration is available for the FOS/RD, where similar errors are known to exist, which 
ultimately limits the wavelength accuracy of our composites to about 0.5-1.0 A (see for 
instance Kerber & Rosa, 2000, for details on FOS calibration). 

Composite spectra cover the entire rest frame region between 1000-3000 A while indi- 
vidual source spectra will often cover only a part of it. We applied the following procedures 



1 Here and below S/N is the signal-to- noise ratio per resolution element, measured in the continuum, 
blueward of CivA1549 
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to compute the composite spectra: 

1. Spectra for a specific source covering the same spectral region were averaged together 
(if they were of comparable S/N). 

2. A single spectrum for each object was created by adding the parts covering different 
wavelength regions, with slight adjustments of the continuum levels if necessary. 

3. All spectra were corrected for Galactic extinction following Schlegel et al. (1998). The 
typical extinction correction Ay was about 0.1 - 0.2 magnitude. No correction for the 
intrinsic extinction was attempted. 

4. The spectra were de-redshifted using the most accurate available measurements (M03) 
or by using the latest accurate values available in the literature whenever a source was 
not part of the M03 sample 2 . It should be noted that redshifts from Veron-Cetty & 
Veron (2001); hereafter VCV01, sometimes differ considerably from our best values 
(Az could be as large as 0.02). 

5. Spectra were normalized by their own continuum flux in the wavelength range 1430- 
1470 A, where no strong lines are observed/expected. 

6. Composite spectra were generated by median-combining the averaged spectra of in- 
dividual AGN. This procedure ensures that only the most typical features will be re- 
vealed, and that some unusual features like narrow absorption (often lost in the noise), 
unusually strong (narrow) lines, noise spikes, etc. will be omitted in the composites. 
A median also works better than an average when the individual spectra do not cover 
the same spectral region (see §4 for details regarding line profiles of the composites). 

The composites have been binned to 0.5 A per resolution element and have a typical 
S/N~50 (see Table 1). The composite spectra are presented in Figure 1 with an enlargement 
of the 1100-2000 A region of greatest interest shown in Figure 2. The latter show details of 
the strongest lines such as Lya, CivA1549, and Cin]A1909. The lines that we identified and 
measured (§3.1) are labeled. 

We defined spectral types on the basis of source parameter properties in the optical 
plane of El (Paper I). This approach differs from the one employed in analysis of the SDSS 
quasar sample (Richards et al. 2002) where sources were binned a posteriori in terms of the 



2 M03 did not publish most of their refined redshift measures based on [OIII]A5007 or the narrow peak of 
H/3 (see Zamanov et al. 2002); they are available from the authors on request. 
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velocity shift between CivA1549 and MgllA2800. We are able to create useful composites 
for five El bins with fixed FWHM(H/9bc) an d RFeii ranges, exactly as done in Paper I. The 
values of FWHM(H/3bc) and R.Feii that provide an optical bin assignment were taken from 
M03. In order to maximize the number of sources in each bin we also include additional 
spectra for sources (~ 20% of the total), not part of the M03 sample. This supplemental data 
involves: 1) sources with reasonable quality archival CivA1549 data and 2) H/3 spectra from 
literature sources with high enough S/N to permit assignment into one of the predefined El 
bins. Details on bin assignation for individual sources that are not part of the M03 sample 
will be provided in a future work. 

Table 1 provides parameter measures and statistical properties for the spectral bins. 
The format is: Column 1 - spectral type designation following Paper I; Col. 2 - range of 
FWHM(H/3bc); Col. 3 - range of RfcIi; Col. 4 - number of sources used for the composite 
spectrum (number of RL source in parenthesis). Each bin typically involves 15 - 30 sources 
resulting in a S/N in the continuum typically of 40 - 60 (Col. 5). Col. 6 and 7 present the 
average absolute V-magnitude (computed from VCV01; Hq = 50 km s _1 Mpc" 1 ; go — is 
adopted) and redshift for sources used in the bin composite followed by the sample standard 
deviations. We also present the average equivalent widths (EW) of [Olll]A5007 and H(3 for 
the objects of the sample - Col. 8 and 9 respectively. 

There is a clear trend between < My > and z, in the sense that Pop. A sources show 
lower mean redshift and luminosity. This is a known bias in the M03 dataset. M03 is a step 
towards our goal to develop a sample complete to my=16.5. At this time the radio-loud 
part of our sample is apparently complete (the standard V/V mSuX completeness test (e.g., 
Peterson 1997) gives ~0.8 compared to 0.3 for RQ sources). The radio loud sources tend 
to be more luminous (due in part to relativistic beaming) resulting in selecting more high 
redshift sources. The El optical parameters show no dependence on luminosity to at least 
z < 1 and m v < 16.5 (Sulentic et al. 2004). We will consider later possible UV line vs. 
source luminosity correlations (see §3.3 and §4.3). 

We also computed the sample standard deviation for the spectra in each bin. Deviations 
from the composites at la confidence level are found to be about 20% and show no specific 
features after being normalized by the composites themselves. 
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3. UV Lines in the El Context 

3.1. Felluv subtraction and line measurements 

Table 2 presents rest-frame EW and FWHM measures for UV lines in the composite 
spectra. El bin designations are listed at the top of the Table. Uncertainty estimates are 
given in parenthesis following each measure. EW uncertainties are about 10-20% for the 
strongest lines (W>10 A) and approach 50% for the weaker lines. FWHM uncertainties are 
generally below 20-30% with the level of uncertainty depending on S/N of the composites 
as well as line blending (e.g. Si IVA1398) in some cases. All measurements were carried out 
after Felluv and continuum subtraction except that these corrections were not applied to 
MgllA2800. 

The Felluv template covering the 1000 - 2000 A region (where the most interesting 
lines are) was created with spectra for I Zw 1 (see Marziani et al. 1996, M03 for details). 
The template was intensity scaled and broadened in order to best match the observed Fell 
pattern in the composite spectra. The results are summarized in the last row of Table 2. 
Figure 3 illustrates the deblending results for two particularly important groups of lines: 
Lya + NVA1240 and Al IIIA1860 + Si III]A1892 + Cm]A1909. The first blend is important 
because NvA1240 is the only high-ionization line that shows no significant Baldwin effect 
(Sulentic et al. 2000a, see also §4.5). The importance of the second blend involves the fact 
that the intensity ratio Si III] A1892/ Clll]A1909 may be sensitive to density (n e ) and/or 
may be affected by metallicity trends along El (M01; §4.3). We carried out a deblending 
procedure by assuming a Lorentzian profile for both lines in the blend and allowing for 
different component intensities and FWHM. A possible contamination by Felll blends may 
contribute to some extend to the relatively poor fit of these lines in A3 bin. 

We consider several interesting results that can be summarized as follows: 

• OviA1034, Lyo, CivA1549, Clll]A1909 and MgllA2800 show a systematic increase in 
equivalent width as one proceeds from extreme population A (bin A3) towards the B 
bins. The most dramatic change is actually seen in MgllA2800, CivA1549 and also 
OviA1034, where bin averaged equivalent widths increase by almost a factor of four (a 
factor ~10 is seen in CivA1549 if individual sources are considered). While we adopted 
CIV centroid shift at FWHM as an El parameter other CIV measures follow the El 
sequence as well. 

• Contrary to the above lines NvA1240 shows almost constant equivalent width (f»20 - 
25 A) across the El sequence. Other lines show no clear trends within measurement 
uncertainties. 
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• The Si III] A1892/ Clll]A1909 intensity ratio decreases monotonically by a factor of four 
as one proceeds from bin A3 to Bl + (from 0.8 to 0.2). See §4 for possible explanations. 

• We find a clear positive trend between FWHM and ionization potential (x) for the 
principal broad lines. The correlation coefficients are in the range 0.6-0.8. FWHM 
typically scales linearly with x U P to x ^ 50 eV and saturates beyond that value. The 
scatter in this relation appears to be smaller for A-bins. B-bins suggest that some LIL 
such as MgllA2800 and Lja deviate from the relationship with larger than expected 
values of FWHM. The linear part of the relation can be roughly described by the 
following expression: FWHM — Vq + ax, where a ~ 60 - 80 and Vo from 1000 to 
3000 km s" 1 from bin A3 to Bl + . The saturation occurs at about 4000 km s" 1 for 
A-bins and at about 5000 - 6000 km s" 1 for B-bins. 



3.2. CivA1549 Profile Decomposition 

The CivA1549 line deserves special attention as the best studied, and easiest to study, 
HIL. The CivA1549 intensity and relative lack of contamination by satellite lines, motivated 
its selection as the representative HIL in the El parameter space. Profile centroid shift 
at FWHM intensity was chosen as one of the most important El parameters (Sulentic et 
al. 2000b). The blueshifted component of the CivA1549 line seems to disappear in sources 
with FWHM H/5bc >4000 km s _1 providing further support for the pop. A-B distinction. 
This fact and other previous works make it clear that both HIL and LIL cannot be treated 
as single component lines. A large part of the confusion in past broad line studies stems 
from this fact. The simplest possible decomposition of our composite CivA1549 profiles 
was attempted after the following processing steps. 1) The continuum was subtracted by 
fitting a low-order Lagrange polynomial over the region 1000-2000A. 2) A broadened and 
scaled Felluv template was subtracted (§3.1, Table 2). Felluv emission is not strong in any 
of the bin composites however a reasonable correspondence between the optical and UV 
Fell emission strengths can be traced in the sense that bin A3 shows the strongest Felluv 
emission. 3) Contaminating lines near CivA1549: Niv]A1486 and blended He IIA1640 + 
Olll]A1663 were fitted with Lorentzians and subtracted from the CivA1549 wings. 4) A 
narrow component was subtracted from the CivA1549 profile (see below). 

Results of the CivA1549 decomposition are given in Figure 4 and Table 3. Results from 
modelling our H[3bc composites (Paper I) are useful here. We identified three reasonably 
unambiguous H/5 components in the context of El and the Population A-B hypothesis. 
These identifications, and the adopted profile fits, are largely empirical at this stage. We 
focus on bins A2/A3 and Bl as representative of the two AGN populations or El parameter 
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space extrema. Bin Al is likely a mix of both population A and B sources (or a transition 
between the two). 

• i) NLR 

H/5: We often see an inflection between the BLR and NLR components of LIL H/3 
making profile decomposition less ambiguous (Marziani et al, 1996). Bins A2/A3 show 
evidence for weak H/5nc emission (EW ~ 3 A) in residuals after modelling the broad 
line components (Paper I). A more significant NLR component is found in the majority 
of Pop B sources < EW >~ 5 — 6 A with a large scatter). If we assume a fixed ratio 
for H/3 NC /[Olll]A5007, then the Table 1 averages for EWQOm] A5007) imply a very 
weak H/9nc contribution in bins A2/A3. This ratio is about 0.11 for pop. B sources 
which would imply EW(H/9nc) — 1 A for bins A2/A3 which is lower than what we 
actually subtracted except in the "blue outlier" sources where [Olll]A5007 is extremely 
weak and blueshifted (Zamanov et al. 2002). 

CivA1549: In our opinion it is clear that a narrow CivA1549nc component exists in 
many, especially Pop. B, sources because we can see a profile inflection (Sulentic & 
Marziani 1999). The situation for CivA1549 is more ambiguous than for H/9 B c because 
of the lower intrinsic resolution (and often S/N) of the CivA1549 spectra which less 
often show the inflection that can guide a profile decomposition. Either Gaussian or 
Lorentzian fits to the broad CivA1549 profiles of bin Al, Bl and Bl + composites show 
a narrow residual. Although one cannot be sure that it is from an NLR component, 
that is the simplest interpretation and is empirically justified by many cases where 
NLR emission is unambiguously observed. New evidence for an NLR component in, at 
least, pop B RL quasars is given in §4.4. We therefore argue that an uncertain NLR 
subtraction is better than no subtraction for the Al and B-bins. No NLR component 
was subtracted from bins A2/A3. The NLR profile width was assumed FWHM< 2000 
km s _1 (taking into account the fact that CivA1549 is a doublet). Photoionization 
models suggest EW(CivA1549 N c)— 0.1 EW([Om] A5007) as a reasonable estimate (i.e. 
about 2 - 4 A, using the EW([0lll] A5007) values from Table 1 and Paper I). We should 
point out that this EW ratio is rather similar to the flux ratio of the lines since the 
continuum shape does not seem to change significantly along the El bins (Fig. 1). 
Without inter-calibration of optical and UV fluxes, as well as knowledge of the NLR 
extinction properties (e.g. Netzer & Laor 1993), such numbers are only indicative. See 
Sulentic & Marziani (1999) for more quantitative estimates of the effect of a CivA1549 
NLR component on measures of FWHM and line shift. Our adopted estimates of 
EW(CivA1549nc)/EW(CivA1549 B c) show a large scatter (0.0 - 0.5) among individual 
sources in Al and B-bins. 
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• ii) BLR 

H/5: This is the defining broad-line component for Type 1 AGN. The Population A-B 
distinction is currently defined in terms of sources with FWHM B.(3 BC smaller or larger 
than 4000 km s -1 . Pop. A2/A3 bin composites are best fit by an unshifted Lorentzian 
function (Veron-Cetty et al. 2001) while for Pop B bins a Gaussian model works 
better. While the bin composite profiles are unshifted, individual sources can show 
red or blue shifts of up to 4000 km s _1 . Profile shifts larger than about 600 km s _1 
are seen only in Pop. B sources (Sulentic et al. 2000a,b). Double peaked profiles are 
also sometimes (<5%) seen in bins Bl + and Bl ++ 3 . Measures of FWHM H/3 B c can 
be strongly influenced by the NLR component (e.g. L/L Edd and BH mass over- and 
under-estimated, respectively). 

CivA1549: A strongly blueshifted line or line component is seen only in Popula- 
tion A sources. The red half of the bin A2/A3 composite profiles can be reason- 
ably well fit with an unshifted Lorentzian model whose strength is inferred from the 
amount of flux redward of the adopted source rest frame. In this context extreme A3 
sources should have only a weak unshifted component with the bulk of the line flux 
strongly blueshifted. We don't know if this decomposition has any physical mean- 
ing, i.e. whether we have two distinct line emitting regions or a single region, pro- 
ducing a complex profile. There are two admittedly empirical reasons to favor the 
two-component model: 

1. H/3bc is best fit with an unshifted Lorentzian and there is no a priori physical 
reason why some HIL emission could not be produced in the clouds responsible for 
the LIL emission. Current models do not rule out the possibility of significant CIV 
emission from the LIL clouds (Korista et al., 1997), although the accurate quantitative 
predictions require extensive knowledge of the physical conditions there. The derived 
width of the best fit CivA1549 Lorentzian components are in good agreement with 
those estimated for H/3bc (about 2500 km s" 1 , see Paper I). 

2. The fits in Figure 4 are consistent with a systematically increasing unshifted compo- 
nent as one proceeds from Bin A3 to Bl. The blueshifted (bulk of the emission in bins 
A2/A3) would be strongly blue asymmetric and would truncate rather sharply near 
the rest frame velocity, an interesting challenge for physical models. This truncation is 
actually seen in sources like I Zwl and Ton 28 where a single component model is more 
appropriate (Sulentic et al. 2000b, 2001). We cannot distinguish between a Gaussian 
or Lorentzian model for the Pop. B BLR. The choice will affect our estimate of the 



3 Bin B1++ was defined in Paper I for sources with R FoII <0.5 and 8000 km s" 1 < FWHM(H/3 BC )<12000 
km s _1 
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strength of the underlying VBLR. Results from H/3bc fits lead us to favor, by analogy, 
a two Gaussian model for CivA1549 Pop. A sources involving unshifted BLR plus 
redshifted VBLR components. At this point we cannot place much physical meaning 
on these Lorentzian vs. Gaussian fits, they simply reflect the relative sharpness of the 
line peak and the extent of the profile wings. 

• iii) VBLR 

B.(3: There is no consensus whether or not the VBLR is a distinct component (e.g. 
Corbin 1997; Snedden & Gaskel, 2004; Korista & Goad 2004). Only Pop B sources 
appear to show this "red shelf or component that we model as a separate VBLR (see 
references in Sulentic et al. 2000a). We prefer to speak of the VBLR as a distinct 
component for at least two reasons: 1) we often see an inflection between the red wing 
of the best fit BLR Gaussian and the red feature and 2) we have identified sources that 
only show a VBLR component (Sulentic et al. 2000c). Paper I composites suggest 
that the best current description of the components involve: i) a BLR with zero mean 
shift and FWHM ~ 5000 - 6000 km s _1 as well as 2) a VBLR with mean redshift of 
about 5000km s" 1 and FWHM~10000 km s -1 . 

CivA1549: No evidence for a significant VBLR component is seen in population A 
bins paralleling the results for H/3bc- A double Gaussian decomposition is consistent 
with both BLR and VBLR CivA1549 components in Pop B bins 4 . FWHM of the 
components are in rough agreement with those for LL3. Interestingly CivA1549 VBLR 
line shifts appear to be lower (~ 1000 km s _1 ) than those modeled in the H/3bc 
composites, however this may be a result of the fact that the two-component fitting is 
not necessarily a unique operation. 

Furthermore line shift values derived for such an ultra-broad component must be con- 
sidered very uncertain especially given the low S/N of many HST spectra. 

While we currently separate population A and B sources at FWHM(H/3bc) — 4000 
km s" 1 , if the distinction has any meaning then it will likely involve a complex combination 
of parameters. There are two possible interpretations of the Pop. A-B concept: i) two 
distinct AGN populations with different central source geometry and kinematics involving 
perhaps a critical Eddington ratio or ii) a single "main sequence" involving the majority 
of AGN and driven by a gradual change in the ionization parameter/Eddington ratio (see 



4 Although we fit Al profile with a single Lorentzian, we cannot rule out the possibility that (some) Al 
sources have both - a weak VBLR and a weak blue-shifted (wind) component, which along with a central 
unshifted component mimic a single Lorentzian. 
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M01). The "elbow" shape of the source occupation in the optical plane of El is an argument 
in favor of the former interpretation (i.e. it is not a monotonic correlation). In any case the 
shape suggests that bin Al will likely blend the parameter space properties of Pop. A and B 
sources. Figure 4 and Tables 1, 2 support that view. A few sources in bin Al (e.g. NGC 1566 
and Mrk 493) show very strong NLR emission while others little or none. Our CivA1549 
profile analysis supports the Population A-B idea in the sense that we find CivA1549 line 
components unique to each (Pop. A - blueshifted BLR, and Pop B - redshifted VBLR) but 
mutually exclusive. 

Our results appear to be robust for sources with z<0.8. Richards et al. (2002) created 
similar CivA1549 composite spectra using part of the much larger and higher mean redshift 
(< z >~ 1.8) SDSS quasar sample. Their spectral resolution is similar to our own but the 
S/N of individual spectra are usually much lower. They binned spectra showing CivA1549 
into four groups (each involving almost 200 sources) based on the centroid shift of CivA1549. 
We lack LIL data that would allow us to assign unambiguous population A and B tags to 
these bins but the two extreme groups (Richards et al. composite D: most blueshifted 
and composite A: least shifted) should, more or less, correspond to our A2/A3 and B1/B2 
samples respectively. If this statement is true then El predicts SDSS composite A sample 
should contain many radio-loud sources and composite D very few. We also note that 
SDSS composite A is more sharply peaked than composite D, consistent with a stronger 
NLR component. We applied our profile decomposition procedure to the SDSS A and D 
composites kindly provided by G. Richards and the results quantitatively verify the previous 
statement about their similarity to our lower redshift composites. The profiles are also 
consistent with the presence of a growing (from D to A) unshifted BLR component that 
further supports our two component fits to the bin A2/A3 composites. Richards et al. 
(2002) interpret the overall shift of CivA1549 as the result of "missing" flux on the red side 
of the profile. We think that the body of evidence, within the El context, does not support 
this interpretation, e.g. the blue wing of CIVA1549 in extreme pop. A sources does not 
resemble the blue wing of CIVA1549 in the broader pop. B sources. Similarly the blue 
wing on SDSS composite D is steeper than that for composite A. If we relate missing flux 
to obscuration then it would be surprising to find the strongest soft X-ray emission from 
sources with the most extreme blueshifted CivA1549 profiles since both should hardly be 
related. 

SDSS composites cannot be used to help us decide whether there is a continuous "main 
sequence" in El or two distinct populations. Caution is needed in interpreting the SDSS 
binning because we do not know what the SDSS bin definition means in an El context, 
especially when MgllA2800 (wings) may be strongly affected by Fell emission. If MgllA2800 
follows the line shift behavior of H/3bc large relative shifts CivA1549 - MgllA2800 can occur 
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in pop. A because of CivA1549 blueshifts and in pop. B because of possible MgllA2800 
redshifts. If CivA1549 blueshifts (relative to rest frame and not MgllA2800) remain a pop- 
ulation A phenomenon at high redshift then the SDSS results suggest that the fraction of 
population A sources is larger at high redshift (75% vs. our estimate that 60% of quasars 
are population A; Marziani et al. 2003b). This would be consistent with our suggestion that 
Pop. A sources represent the "seed" population of young quasars in the Universe (Sulentic 
et al. 2000a; Mathur, 2000; Constantin & Shields, 2003). 

3.3. RL — RQ Differences and Luminosity Effects 

We can also compare RL vs. RQ properties using suitably constructed composites. 
Selection of RQ and RL subsamples from a single bin will ensure that we compare quasars 
with similar physical and observational properties. An extreme comparison of Pop. A vs. 
B will not be possible using a single bin. We would have to compare an extreme pop. A 
bin with an extreme pop. B one. Instead this will be a comparison between RL sources 
and RQ sources with optical spectra indistinguishable from them. This comparison offers 
hope to isolate subtle effects that might be direct manifestations of radio loud activity. We 
must use a B bin for this comparison because that is where most RL sources are found. We 
choose Bl which is the most densely populated and where RL and RQ objects are almost 
evenly represented. Radio (6cm) and optical V-band fluxes were derived from data in the 
VCV01 catalog. A radio-index (Kellermann et al. 1989) Rk was computed and all sources 
with log(i?K)>2 are considered radio- loud (see Sulentic et al. 2003). Figure 5 shows the 
RL and RQ composites, as well as RL-RQ residuals, for both Lya and CivA1549. We see 
evidence for differences: 1) most significantly, stronger (EW~5A), NLR emission in the RL 
sources and 2) a possible excess redshifted BLR (or VBLR) emission in RL sources (this is 
also present in Figure 7 A-D composite difference of Richards et al. 2002). The former is 
not a surprise (see e.g. Xu et al. 1999) while the latter, if real, may be an important clue. 
Since we know that some bin Bl RQ sources also show NLR emission result 1) reflects a 
lower limit for the mean NLR contribution in Bin Bl RL sources. A comparison shows also 
that the RL sources from this bin have on average stronger [OIII] and H/Snc components 
(typically by about 30%) than the RQ sources. There is no doubt that RQ sources exist with 
VBLR components as strong as the strongest observed among the radio-louds (e.g. Sulentic 
et al. 2000c) so further tests of the red residual are needed. 

In order to explore any luminosity correlation with composite spectra we again restrict 
ourselves to bin Bl. Sources were separated into brighter and fainter groups using My = —25 
as the reference value. This value separates the bin Bl sample into two approximately even 
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groups. Composite spectra (covering 1000-2000A) and differences are shown in Fig. 6. One 
can see a difference in the broad components of all the strongest high ionization lines ( 
Lya, CivA1549, Si IVA1398 and Cin]A1909) in the sense that the equivalent width is larger 
for the lower luminosity composite spectrum (the Baldwin effect? see §4.5 for a thorough 
discussion). This difference was not seen in the RL vs. RQ comparison. Fainter, lower 
ionization lines do not show any difference, furthermore they are weak enough that any 
difference would likely be undetectable. 



4. Discussion 

4.1. The Relevance of the Median Composites 

The use of composite spectra raises an immediate question. Do median composites yield 
typical quasar spectra for each bin? For example, is a symmetric profile in a median com- 
posite the result of equal numbers of blue and red asymmetric profiles? A visual inspection 
of our sample spectra shows, e.g. for B-bins, about 60% of sources with a more or less sym- 
metric CivA1549 profile, while the remainder show 10, 20 and 10%, respectively red, blue 
and indeterminate (due to low S/N) asymmetries. This suggests that median composites 
represent the most typical quasar spectrum and that the changes between our spectral bins 
are real. If our bin populations had > 100 objects then we would use weighted average com- 
posites, but medians appear to be a reasonable approach at this time. It is also important 
to point out that bin sizes were chosen to include sources that are statistically similar within 
the average parameter measurement errors. 

Another important question is whether or not an ensemble of line profiles with a wide 
distribution of widths, strengths and asymmetries will yield a representative profile of the 
same type after combination. We see that the many of the lines (especially the strongest) are 
better fit by a Lorentzian whatever that means physically. Is it possible, for instance, that 
an ensemble of Gaussian profiles will generate a Lorentzian-like shape in composite spectra 
or the converse? This might lead to incorrect interpretation, and hence modeling, of line- 
forming regions. In order to test this possibility we ran a simulation by combining together 
about 100 artificial line profiles with a wide (normal) distribution of widths and intensities, 
resembling the real situation to some extent. We performed simulation with both Gaussian 
and Lorentzian profiles, creating median and average composites from the individual pseudo- 
spectra. The results suggest that profile type is preserved, i.e. a composite of Gaussians or 
Lorentzians will preserve their respective profile forms in a composite spectrum. 
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4.2. Inferences Prom CIV Profile Decomposition 

Proper decomposition of emission lines in composite spectra along the El sequence 
offers insights into changes in BLR geometry and kinematics. Several results from this and 
previous papers provide interesting clues. 

Asymmetric blueshifted HIL emission is seen almost exclusively in population A sources. 
In the most extreme blueshift is also observed in the unusually weak [Olll]A5007 

features (Zamanov et al. 2002; Marziani et al. 2003c). One simple interpretation that 
immediately suggests itself involves some kind of outflow from an optically thick emitting 
region into an optically thin NLR. In a scenario where the broad lines arise from an accretion 
disk this could involve a high ionization disk outflow or wind. It would be strongest in 
population A sources and would gradually weaken or even disappear as one proceeds along 
the El sequence toward population B. Along the same sequence an unshifted/symmetric 
Lorentzian HIL component (which may be analogous to the dominant LIL component) may 
be growing, in exactly the opposite sense of the blueshifted feature. In fact some theoretical 
models predict Lorentz-like profile wings (e.g. Penston et al., 1990; Dumont & Collin- Souffrin 
1990). The profile may become more Gaussian in the population B region. We have argued 
in a pre-El context that HIL vs. LIL line shifts appeared to be completely uncorrelated for 
RQ sources and marginally correlated for the RL minority (Marziani et al. 1996). Composite 
spectra, binned in an El context, suggest that CIV and H/5bc ma Y show analogous emission 
components except for the blueshifted HIL feature seen only in pop. A that are almost 
always RQ AGN. 

A very broad and redshifted VBLR component is only seen in population B sources and 
may be present in both HIL and LIL. Since NLR-BLR-VBLR inflections are rarely seen, 
there may be less spatial/kinematic decoupling between these regions for CIV than is seen 
for H/Sbc, where clear inflections are common. That is far from proven because we also 
see LIL profiles with little or no inflection between components. Inflections may be driven 
by kinematic and/or orientation. There are some indications that the region producing 
CivA1549nc emission is closer to the center that the region producing Hj3 NC (Sulentic & 
Marziani 1999). The ClvA1549 VBLR may be most simply interpreted as arising at the 
inner edge of the classical BLR in analogy to H/3 B c VBLR (Marziani et al. 1993; Sulentic 
et al. 2000c). The lower redshift of the CivA1549 Pop. B. sources VBLR may also indicate 
that a high ionization wind is still present in the VLBR since the ClvA1549 component is 
still blueshifted with respect to the corresponding B./3 component. 
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4.3. BLR Structure and Kinematics 

So far there is no single picture capable of explaining the overall properties of the broad 
line emitting region. Currently fashionable models involving line emitting clouds near an 
optically thick accretion disk suffer a number of difficulties. If the motions of the clouds are 
Keplerian then they will collide with the disk and will likely be destroyed on short time 
scales. If one adopts an outflow/inflow model then the entire profile will show large red or 
blue-shifts because, the standard optically thick disk will obscure emission from the far side 
of a flow, while in reality large shifts in line profiles are rare. 

Alternatively, emission lines originating from a geometrically thin, optically thick disk 
(e.g. Dumont & Collin- Souffrin 1990) will show extremely small FWHM when observed 
face-on. This situation may be relevant for pop. A sources but latest work suggests that the 
situation is different for RL and, by analogy, other pop. B sources. The most jet-aligned 
superluminal RL sources, which, should have negligible disk rotation contribution to LIL 
FWHM, show FWHM~3000-4000 km s" 1 (Rokaki et al., 2003; Sulentic et al. 2003). It is 
however true that core-dominated (CD) and superluminal sources with % —>■ 0° show sys- 
tematically narrower profiles than other (FRII) RL AGN (Sulentic et al. 2003). CD sources 
are displaced towards lower FWHM(H/9bc) in the optical El plane, which is consistent with 
orientation predictions. However, even these near pole-on sources lie well above many RQ 
sources with similar RFeii, i-e., they have larger FWHM(H / 9 BC ) for the same R Fe ii (Sulentic et 
al. 2003). The results for radio loud sources suggest that we need a considerable velocity dis- 
persion in the vertical direction to account for observed LIL widths (FWHM(H/? B c)~ 3000 
km s _1 ). 

The above results might be telling us that the bulk of line emission in RL and Population 
B RQ sources does not arise from a standard optically thick accretion disk especially when 
we observe FWHM LIL/HIL ~ 3-4000 km s _1 in sources reasonably interpreted as almost 
disk face-on. One candidate for the line emitting region involves the outer self-gravitating 
part of the disk (Collin & Hure, 2001; Bian & Zhao, 2002). The disk is expected to become 
gravitationally unstable at some radius Rsg, where it will break up into rings or discrete 
self-gravitating clouds. The vertical velocity dispersion of such clouds may be amplified by 
cloud-cloud encounters. Some studies suggest that the disk can become self-gravitating at a 
distance from the central black hole as small as 100-200 R G , where R G is the Schwarzschild 
radius (see for instance Collin & Hure, 2001). A detailed theory for a self-gravitating disks 
is not yet available, however it is expected that the radius of the self-gravitating region 
will depend on the central black hole mass and Eddington ratio (and two related and very 
important parameters: viscosity, a, and disk opacity). There is no consensus on the correct 
dependence of R$g on parameters like M B h or Eddington ratio m = L/L Edd (here we assume 
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that the accretion efficiency (77) does not change with luminosity, otherwise m = M/M^dd, 
where M Edd = L Edd /(r]c 2 ). 

A very simple approach for exploring this idea involves application of the Toomre sta- 
bility criterion to a standard Shakura-Sunyaev disk (Shakura & Sunyaev 1973). This will 
give an estimate of the radius where the stability parameter (Q) becomes approximately 
1. Different models of the self-graviting region often predict different Rsg dependence on 
the accretion parameters because of different initial assumptions about the outer regions 
of the disk (i.e. opacity, pressure, etc). Bian & Zhao (2003) obtain in a recent work: 
.Rsgoc m~°' 49 Mg^' 16 based on standard accretion disk solutions assuming dominant Kramers 
opacity. They assume that it is equal to the radius of the BLR /?blr (here and below i?sc is 
expressed in R G ). Collin & Hure (2001) adopt a numerical approach with realistic opacities 
and get a complex R SG dependence roughly approximated by Rsgoc m°' 16 Mgg 46 for the 
region log(M B H)=7 - 9, m= 0.01 — I. Following these authors, for a BH mass M=10 7 M 
and Eddington ratio ~1, .Rsg- 10000 R G , while for M=10 9 M and m~0.01, R S g is less 
than about 500 R G . 

On the other hand we can consider " superluminal" radio- loud sources with the narrowest 
RL line profiles: FWHM(H/3 B c)^ 3000 km s _1 , and the narrowest Pop. A sources, with 
FWHM(H/5bc)^ 1000 km s _1 , also likely to be observed at i ~ 0°. Assuming orbital motion 
with Keplerian angular velocity (Qk), one can write for these apparently face-on sources: 

FWHM ~ Av ~ istt K R SG (1) 

where Av is the vertical velocity dispersion, assumed to be proportional to the Keplerian 
velocity by a factor of v. A reasonable guess for v is about 0.1-0.2. One can easily show 
that FWHM of 3000 km s" 1 implies a radius of about 5000 R G while 1000 km s" 1 implies 
500 i?c, values quite close to the size of the self-graviting region assumed above. Even if the 
agreement is perhaps fortuitous, and the exact Rsg dependence is far from certain, it is rea- 
sonable to conclude that R SG (~ -Rblr) may be smaller in Pop B sources by a factor ~ 10, 
and that this may leave a very "small" emitting surface for any standard optically-thick 
geometrically-thin disk component. Part of the line profile may be produced in the disk if it 
is illuminated by a geometrically thick inner region or by disk warping (e.g. Bachev, 1999) 
and if it is not extremely hot (which would completely ionize the lower-ionization species). 
On the other hand, the presence of a blue-shifted component can be most simply explained 
as the signature of an outflow from an optically thick disk. The disk-wind scenario can pro- 
duce blue-shifted high-ionization lines as already suggested in several works (e. g. Murray 
& Chiang 1997; Bottorff et al. 1997). If our interpretation is correct this implies that the 
wind strength increases toward the most extreme population A sources (bins A2/A3 with 
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strongest RFeii)- Population B sources may have a much lower, or may not develop, a wind 
component. Following these considerations we propose a simple BLR scenario (see also Fig. 
1 from Collin & Hure, 2001): 

• A "hot" innermost accretion region, located inside some transition radius i?x r ^ 50 — 
100/?g)- Our data leave the structure of this region open to speculation, although we 
may reasonably assume that the continuum region (or part of a geometrically thick, 
optically thin advection-dominated accretion flow (ADAF; see Narayan, Mahadevan & 
Quataert 1998 for an exhaustive review) can heat the outer line-emitting region. 

• An optically thick, geometrically thin disk whose extent decreases by a factor ~10 
from Pop. A to Pop B. The heated disk produces a wind that is responsible for the 
blueshifted CivA1549 component. The outer region of the disk may emit H(3bc an d 
other LILs as well in Pop. A sources. The disk may become almost negligible in Pop B 
sources because Rsg may become as small as 100 - 200 Rq (e.g. Bachev & Strigachev, 
2004). 

• A self-gravitating and fragmented zone where line emission might arise in dense gas 
clumps. This region may be the entire BLR for Pop. B sources while it may contribute 
mainly to LIL emission in Pop. A objects. 

Within the context of this model + El we can eventually explain several observational 
results from this paper as well as from previous works: 

• The different vertical velocity dispersion in the self-gravitating zone for Pop. B (closer) 
and Pop. A (farther) sources accounts for different line widths at % — > 0°; 

• After subtraction of an unshifted/symmetric Lorentzian component the blueshifted 
residual shows reasonable persistency in shape - the intensity may change but not the 
overall shape. The blue component becomes marginal at bin Al. This behavior is 
consistent with the idea that the blue component comes from a disk-generated wind, 
where the kinematics (velocity field) is strongly affected by the detailed properties 
(e.g. magnetic field structure and terminal wind velocity). The overall shape of a line 
arising in a disk wind will be less dependent on the radial distance where the wind 
originates. In other words, our model provides an easy explanation why population A 
sources tend to show decoupled HIL-LIL components while Pop. B sources show more 
similar CivA1549 and H/? B c profiles. In the most extreme cases, there are sources 
where the entire CivA1549 profile resembles the blue component alone (e.g. I Zw 1, 
Ton 28; Sulentic et al. 2000b). 
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Leighly (2004) presented an interesting explanation of the apparent anticorrelation 
between the strengths of the central and the blue-shifted (wind) components. The wind 
that gives rise to the blue-shifted component might filter the EUV ionizing radiation 
that would otherwise reach the outer, self-graviting parts (where the unshifted line 
component might originates). The wind cannot however filter effectively the hard 
X-rays that are reprocessed in LIL emission in the same region. 

An alternative approach is to ascribe the entire profile of CIV in A2/A3 bins to an 
outflow (wind). In such a case the redward emission is seen either because the disk 
is (partially) transparent to CIV photons or because of higher inclinations, resulting 
in significant projection of the Keplerian/radial component of the wind velocity onto 
the line of sight. Both possibilities seem somewhat unlikely, at least for these extreme 
bins. 

• There are trends between x and FWHM (e.g. Tytler & Fan 1992). Pop. A sources can 
be explained with a split BLR: inner HIL and an outer LIL emitting regions. Since 
we expect that gas motion will retain its circular velocity in the wind and will acquire 
an upward and outward velocity component, the velocity dispersion in the wind will 
be higher than for a Keplerian velocity component in the disk. LIL emission from the 
outer, fragmented region will therefore be narrower. In Pop B sources we expect a 
"normal" BLR stratification: HIL emitted closer to the central continuum source than 
LIL. In general, we expect that both HIL and LIL will be emitted rather close to the 
central continuum source, accounting for the rather high ionization spectrum. This 
outline requires testing with photoionization models. 

Reverberation mapping indicates continuum-response times of Cin]A1909 consistent 
with the one of B.(3 BC (at least in NGC 5548; Clavel et al. 1991). This result suggests that 
Cin]A1909 is primarily produced in the low-ionization emitting part of the BLR. Under the 
assumption of nearly constant [Si/C] abundance ratio, the behavior of the Si III] A1892/Cin] A1909 
intensity ratio along the El sequence can then be understood in terms of increasing electron 
density toward type A3 (M01) in the LIL-emitting BLR. The critical n c of the Cm]A1909 
line is ~ 10 9 ' 5 cm -3 , much lower than that of the Si III]A1892 line, ~ 10 11 cm~ 3 . There- 
fore, for n e > 10 10 cm -3 , the Cin]A1909 line should be more and more collisionally sup- 
pressed with increasing n e , while the Si III] A 1892 line should become favored. We also 
observe a consistent increase in the strength of the Al IIIA1860 line, with the intensity 
ratio Al IIIAI86O/C111] A1909 reaching ~ 0.8 for spectral type A3 (n c ~ 10 11 cm~ 3 ). The 
Al IIIA1860 line is expected to strengthen with increasing electron density (Baldwin et 
al. 1996; Korista et al. 1997). However, as pointed out by the referee, chemical abun- 
dance may not be constant along the El sequence. We do not have information on the 
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[Si/C] and [Al/C] abundance ratios. Intensity ratios involving nitrogen lines suggests that 
super-solar metallicity is appropriate for extreme population A (i.e., A3) sources. From 
Table 2, we have NvA1240/CivA1549^ 1 and NvA1240/He IIA1640^3 for A3, while for 
Al NvA1240/CivA1549^ 0.25, and NvA1240/He IIA1640«1.5. It is intriguing that we are 
observing again, for spectral type A3, line ratios that are more frequently found in high- z 
quasars (Hamann & Ferland 1999). Highly supersolar nitrogen enrichment may have been 
produced by star formation with a top-loaded initial mass function (Hamann & Ferland 
1992; §4.5). It is possible that changes of line ratios from Pop. B to Pop. A and along the 
A-spectral type are due to both an n e -increase and a metallicity increase. Quantifying these 
effects is a challenge for further work. 



4.4. Radio Loud Sources: Black Hole Spin and Host Morphology? 

Comparison of RL and RQ line profiles (§3.2) in the El context should begin with a 
comparison of population A and B sources. The differences that we find between Pop. A 
and B HIL CIV composites are even larger than for LIL B.(3 B c- We can say with some 
confidence that sources showing low EW, blueshifted and blue asymmetric CivA1549 have 
low probability (P<1%) to be radio-loud. Pop.B sources have a reasonably high probability 
to be RL (P~50%; Marziani et al. 2003b). Since not all pop. B sources are RL we 
focus (§3.3) on a particular bin (Bl) within population B where we can compare CivA1549 
properties for RQ and RL AGN that appear to be spectroscopically indistinguishable. Our 
RL-RQ difference spectrum is shown in Figure 5. We see evidence for a significantly stronger 
NLR component in pop. B RL sources. This might be explained if a significant part of the 
narrow line emission is produced in interactions between a radio jet and the ambient medium, 
forming presumably a narrow-line-bright cocoon around the jet (e. g., Steffen et al. 1997; 
Axon et al. 1998). Aside from a possible excess of redshifted broad line gas in the RL sources, 
Pop. A RL and RQ line profiles look very similar as was also found for H/3bc (Marziani et 
al. 2003b). The VBLR component, which might be related to this redshifted gas, is seen in 
both RL and RQ Pop. B AGN (Sulentic et al. 2000c) and, therefore, is not a RL signature. 
One is tempted to consider host galaxy morphology as a major difference between the RL 
and RQ sources in bin Bl. RL sources tend to be hosted in ellipticals (Pagani et al. 2003) 
while a typical pop. B RQ source involves a classical Seyfert 1 like NGC 5548, hosted in a 
spiral. Perhaps the Pop. B emission line spectrum indicates that RL activity is possible as 
far as physical conditions near the BH are concerned but host morphology dictates whether 
it happens or not. If it were that simple however we might expect to see many frustrated RL 
sources involving trapped core emission (unless it is self quenching under that circumstance). 
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The bin Bl results emphasize the fact that Pop. B RQ and RL sources are spectro- 
scopically indistinguishable as far as broad lines are concerned. This may be telling us that 
radio-loudness is not directly dependent on, or strongly influencing, BLR structure (and 
therefore - the structure of the accretion disk, at least at r > 100 Rq)- This makes a cir- 
cumstantial case for the role of black hole spin, which is unlikely to have any influence on 
the accretion disk structure outside 10-100 Rq where the UV line emission likely originates, 
at least if the Bardeen-Petterson effect (Bardeen & Petterson, 1975) is not in play. The 
black hole spin can, however, play a significant role in jet production (e.g. Blandford & 
Znajek, 1977). The fact that the majority of RL sources occupy B-bins may be associated 
with a change of the innermost accretion disk structure (from a thin disk to an ADAF?, 
Rozanska & Czerny, 2000, and references therein) that is probably triggered by the lower 
accretion rate (or by some function of the black hole mass and the accretion rate). If true, 
the accretion disk structure change (from thin disk to ADAF) will be a necessary but not 
sufficient condition for radio-loudness (i.e. the ultimate trigger will be the black hole spin, 
subject to the environmental caveat mentioned above). In this scenario only B-bin objects 
which contain a rapidly spinning black hole would be able to produce powerful jets. 

4.5. A Local "Baldwin Effect" Driven by the Eddington Ratio? 5 

Measurement of broad line EW's reveals both a trend and the lack of a trend (see also 
§3.1): 

• Lya, CivA1549 MgllA2800 and OviA1034 show a systematic EW decrease by factors 
of 2 - 4 as one proceeds from bin Bl + to A3. Table 2 shows the same trend for other 
lines but caution is needed when interpreting blended lines and data for lines with 
EW<5A that have estimated uncertainties >30% at the 1-cr confidence level. 

• The NvA1240 line shows no trend with EW(NvA1240) « 20-25 A in all spectral bins. 

These results essentially describe the so-called "Baldwin effect" (hereafter BE; originally 
defined as an anti-correlation between rest frame EW CivA1549 and source luminosity at 
1450A: Baldwin 1977). It is well known that NvA1240 does not follow the BE trend (see 
Sulentic et al. 2000a for references before 1999; Dietrich et al. 2002; Green et al. 2001; 



5 We call the L/LEdd-dependent Baldwin effect also a "local" Baldwin effect since it is defined through 
low-z, low to moderate luminosity quasars and Seyfert 1 nuclei, in opposition to the canonical Baldwin effect 
which emerges considering high-z, high-luminosity quasars 
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Croom et al. 2002). While the original BE showed a tight correlation involving a few tens 
of sources, more recent work (Dietrich et al. 2002; Warner et al. 2003) show it to be a very 
loose correlation that only becomes significant over a wide range in continuum luminosity 
log(AL>) =42-48 ergs s _1 (e.g. Kinney et al. 1990). Considering the large dispersion it is not 
surprising that several studies using small samples found no evidence for a BE (see Sulentic 
et al. 2000a for numerical simulations and a thorough discussion of the issue). 

The BE-like trend found between our spectral bins does not represent a dependence on 
luminosity. The El sequence is independent of source luminosity (Sulentic et al. 2000b) and, 
is more likely governed by the Eddington ratio (M01; Marziani et al. 2003b). In fact Pop. B 
spectral bins contain more high luminosity sources than Pop. A because RL in our sample 
are overluminous, and overrepresented, relative to RQ sources (see also Boroson & Green 
1992). Our BE trend is also not BH mass-dependent (e.g., Warner et al. 2003): the largest 
BH masses are again found in population B (M03). If we compare the average m values for 
our bins (from M03b) with EW(CivA1549) measures from our composite spectra we obtain 
the following approximate relation: log EW(CivA1549) ~ 1.35 — 0.67 log excluding bin 
Bl + . Our models predict that Bl + sources are simply more inclined Bl sources (M01). 
Including Bl + yields log EW(CivA1549)~ 3.22 — 0.39 log -j^ using L/M bin averages from 
M03b or log EW(CivA1549)~ 3.54 — 0.47 log if bin averages are taken from this paper 
(L/M is measured in Solar units), f we consider all sources (RQ and RL) a robust fitting 
method yields log EW(CivA1549)~ 3.23 - 0.36 log ^, which becomes log EW(CivA1549) 
~ 3.40 — 0.41 log for RQ sources alone. 

Parameter dispersion in the Baldwin effect (from EW(CivA1549) ~ 20 - 200A) is well 
reproduced by sources spanning the range ~ 0.02 - 1 in m regardless of the fitting solution 
details. Highest m sources show the lowest EW(CivA1549) and it is a well-established fact 
that NLSyls show the lowest EW(CivA1549) (Marziani et al. 1996; Rogriguez-Pascual et 
al. 1997; Sulentic et al. 2000a) and are likely to radiate close to the Eddington limit. Many 
known NLSyls are relatively low redshift /luminosity sources that tend to blur standard BE 
correlations at the low luminosity end (see e.g., Brotherton & Francis 1999) because they 
show EW(CivA1549) similar to the high redshift quasars (Sulentic et al. 2000a). Sulentic 
et al. (2000a) suggested that the BE might be caused by preferential (probably driven by 
selection effects and/or intrinsic evolution) detection of large m sources at high redshift. 
In a flux limited sample, where L and z are highly correlated, an "evolutionary" (redshift- 
dependent) Baldwin effect is expected and is indeed found in the Large Bright Quasar Survey 
(Green et al. 2001). Given a sample that is complete in terms of m, the BE may survive with 
even larger scatter or disappear depending on the relative importance of evolutionary and 
selection effects (see Dietrich et al. 2002). Deeper surveys may detect larger EW(CivA1549) 
sources at high redshift: see for example the distribution of EW(CivA1549) at z > 4 in 
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Constantin et al. (2002). Intriguingly, they detect several sources with EW(CivA1549)>50A 
at logL m 46.5 ergs s _1 . 

Interpretation of the BE has been highly controversial with suggestions encompassing 
orientation effects (Netzer 1985), selection effects (Sulentic et al. 2000a, and references 
therein), black hole mass (Wandel 1999; Warner et al. 2003), Eddington ratio (Sulentic et 
al. 2000a, Shang et al. 2003; Baskin & Laor 2004). A continuum softening with increasing 
luminosity (Wandel 1999) is a less model-dependent explanation. Warner et al. (2003) 
suggest that the driving parameter may be the central black hole mass. These authors 
computed the BH mass using a modified Kaspi et al. (2000) relation, employing FWHM 
CivA1549 instead of FWHM H/3 and assuming virial motions. This does not demonstrate 
that mass is driving the BE because it is assumed to be M oc L - 7 , and this assumption makes 
the computation in part circular. Furthermore, we question the validity of estimating BH 
masses using CivA1549 (e.g. Vestergaard 2002 and Warner et al. 2003; 2004). While H/? B c 
shows average broad line profiles that can be argued to arise in virialized clouds, CivA1549 
is blue asymmetric and blue shifted in about 60% of sources. In many cases the shift is 
significant and in some - virtually all of the line flux is displaced to the blue side of the local 
rest frame. It is not at all clear that the clouds (or wind?) producing the HIL emission can 
be assumed to be virialized. Finally failure to correct for NLR emission in the CIV profile 
will contribute to the dispersion in BH mass that might be present. NLR emission tends to 
decrease the dispersion BLR FWHM measures (FWHM measures for pop. B become more 
like Pop A). Figure 5 shows that NLR uncorrected FWHM CIV measures will underestimate 
BH mass and overestimate L/L Edd values for RL sources. 

Our work shows that a CivA1549nc exists in, at least, Pop. B sources and that a 
significant part (if not all) of CivA1549 emission in bins A2/A3 arises in a wind. The 
terminal wind velocity could be only indirectly related to Keplerian (i.e. virialized) motions. 
The overall velocity dispersion is even less likely to be related to virial motion in near face-on 
sources, where the accretion disk likely obscures the receding part of the flow. The CivA1549 
line profile in such cases is almost fully blueshifted with respect to the rest frame. While 
any correlation with mass will be little affected over several orders of magnitude (Warner et 
al. 2003), any more fundamental correlation with rh may be entirely washed out (rh range 
is a factor <50; Woo & Urry 2002). Comparing FWHM H/3 and CIV in different spectral 
bins (Pop. B with and without CivA1549nc subtraction) suggests that errors could be at 
least a factor ~10 for some objects (since, for example, the FWHM ratio CIVA1549H/3 is « 
3 in IZwl; Marziani et al. 1996; 2003b). However, a shallow trend between EW(CivA1549) 
and L/L Edd is noted in a large sample of quasars of redshift < z < 5 (Warner et al. 2004), 
which seems to support our results for z <0.8. 
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Another circumstantial element supporting a role for the accretion rate m in the BE in- 
volves the constancy of EW(NvA1240) and the increasing intensity ratio NvA1240/CivA1549 
along the sequence from Bl + to A3. The large values of EW(NvA1240) for all spectral types 
are suggestive of solar or supersolar metallicity (Turner et al. 2003; Bentz & Osmer 2004). 
It is not unreasonable to suppose that N may have been enhanced by vigorous formation of 
massive stars that burn H via the CNO cycle in the AGN circumnuclear regions. The largest 
enrichment may have occurred in relatively young or rejuvenated quasars radiating at large 
accretion rates as NLSyls are thought to be (several NLSyls hosts may be merging dwarf 
or barred galaxies; Krongold et al. 2001; Crenshaw et al. 2003). 

5. Conclusions 

We computed and analyzed composite UV spectra for 5 different spectral bins along 
the main sequence of the El optical parameter plane. Data were analyzed following Paper 
I, which considered similar composite spectra for the optical region involving H/3, our repre- 
sentative LIL. We measured the FWHM and equivalent width of a number of emission lines 
focussing on CivA1549, our representative HIL. Our main results include: 

• We decompose HIL CivA1549 into possible NLR, BLR and VBLR components many 
of which have LIL analogs. We consider them to be distinct line emitting regions. Type 
2 AGN show only NLR LIL and HIL emission. Type 1 AGN are know with pure BLR, 
pure VBLR and, most often, composite BLR+VBLR (with and without NLR) spectra 
if they belong to Population B. A strong asymmetric blueshifted BLR HIL component 
is unique to population A (bins A2/A3) sources. A reasonably strong and redshifted 
VBLR is unique to HIL and LIL in population B (bins Bl and Bl + ) sources. 

• We infer a n e trend along the El sequence using the ratio of the semi-forbidden lines 
Si III] A1892 and Clll]A1909. Under the assumption of a nearly constant chemical 
abundance along El, we find that density increases by about an order of magnitude 
between bins Bl + and A3. However, there are indications that both - density and 
metallicity may increase from extreme B-bins to extreme A-bins. 

• We compare sub-samples of RL and RQ sources in bin Bl. We find an excess of NLR 
and redshifted broad line emission in the RL subsample. The former represents virtual 
proof that NLR line emission is present in the spectra of many type 1 AGN. 

• We advance the following empirically driven model for the structure of the accretion 
flow: 
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1. An inner hot optically thin (possibly ADAF) disk, located inside Rtt < 100 Rq- 

2. A standard thin disk, between and -Rsg, producing most of the optical contin- 
uum. This is the place where a wind may originate and produce the blue-shifted 
CIVA1549 component. 

3. An outer self-gravitating region at r > -Rsg? producing most of the broad-line 
emission. 

i?Tr and Rsg are functions of the accretion parameters. A n eventual discontinuity 
in their behavior as functions of M and m may account for the existence of different 
quasar populations (e.g. Pop. A, B). 

• We suggest that the Eddington ratio is the primary driver of the Baldwin effect rather 
than source luminosity. This suggestion is motivated by the robust BE that we find 
in El which is luminosity independent. This explains why high redshift sources might 
preferentially show low EW and blueshifted CivA1549 if Pop. A sources are the "seed" 
(i.e. youngest high accreting lowest mass BH) AGN. This interpretation is supported 
by comparison of our low redshift sample with first results from higher redshift SDSS 
quasars. The latter show an increase in the fraction of sources with a blueshifted 
CivA1549 HIL - a defining characteristic of Pop. A AGN. 
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Table 1. Definition of spectral types 



Bin 


FWHM(H/3 BC ) 


Rpcii 


Ntot (Nrl) 


S/N 


<M V > (cr) 


<z> (cr) 


EW([OIII]) 


EW(H/3) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


A3 


<4000 


>1.0 


10 (2) 


39 


-23.1 (2.1) 


0.20 (0.20) 


20 (5)* 


50 (54) 


A2 


<4000 


0.5 - 1.0 


14 (1) 


35 


-24.4 (1.6) 


0.26 (0.22) 


18 (9) 


67 (63) 


Al 


<4000 


<0.5 


40 (12) 


74 


-24.3 (2.2) 


0.26 (0.24) 


44 (25) 


95 (89) 


Bl 


4000 - 8000 


<0.5 


54 (28) 


87 


-25.1 (1.7) 


0.36 (0.24) 


46 (26) 


98 (95) 


B1+ 


>8000 


<0.5 


21 (17) 


52 


-25.4 (1.4) 


0.44 (0.22) 


60 (34) 


99 (98) 



*We provide the median values (in brackets) for the last two columns as well. Due to significantly dispersed EW 
distributions and the small number of objects in some bins, the medians may differ significantly from the averages, 
especially for [OIII]. 



Table 2. Line measurements 



LjIiig Id. 




A 3 




A2 




Al 




Bl 




B1+ 




TTW (rr\ 
H> VV \ U J 


r vv xiivi ycr j 




r vv nivi \cr ) 


i2j vv yd ) 


r vv xiivi ^cr ) 


Hi VV \<T ) 


r vv nivi \cr ) 


I2j VV y(7 j 


r vv nivi yu j 


OVI 1034 










17 (4) 


3770 (500) 


21 (6) 


4350 (600) 


18 (4) 


5510 (700) 


CIII 1176 






4 (3) 


3600 (2000) 


1 (1) 


5100 (2000) 


4 (3) 


5350 (2500) 






Lyal216 


48 (15) 


2960 (700) 


92 (21) 


2960 (500) 


101(25) 


3450 (500) 


116(21) 


3700 (400) 


93 (18) 


5420 (400) 


NV 1240 


26 (9) 


3870 (700) 


21 (7) 


4110 (700) 


17 (6) 


3380 (600) 


28 (9) 


5080 (800) 


25 (6) 


7250 (900) 


Sill 1264 


3 (2) 


2250 (1000) 


2 (1) 


1420 (800) 














Sill 1306 


6 (2) 


3530 (500) 


5 (1) 


2290 (300) 


3 (1) 


2290 (300) 


4 (1) 


3670 (500) 


3 (1) 


4820 (1000) 


CII 1335 


2(1) 


1970 (500) 


3 (1) 


1680 (300) 


1 (1) 


2240 (300) 


1 (1) 


2240 (600) 


1 (1) 


3590 (1500) 


SilV 1398 


16 (3) 


4500 (400) 


16 (3) 


3860 (400) 


13 (3) 


4070 (400) 


11 (2) 


3640 (400) 


12 (2) 


6000 (400) 


NIV] 1486 










2 (1) 


2610 (800) 


1 (1) 


3020 (800) 






CIV 1549 


23 (4) 


4160 (400) 


35 (9) 


3290 (400) 


69 (16) 


3480 (400) 


83 (16) 


4260 (400) 


79 (10) 


6190 (400) 


Hell 1640 


8(2) 


5850 (700) 


10 (2) 


4750 (500) 


11 (2) 


4570 (400) 


11 (2) 


4930 (500) 


4(1) 


4750 (600) 


OIII] 1663 


2 (1) 


2880 (800) 


5 (2) 


4140 (700) 


4(1) 


3060 (500) 


7(2) 


4320 (700) 


5 (2) 


3600 (700) 


NIII] 1750 


4(3) 


2050 (1000) 










4 (3) 


4410 (1500) 






A1III 1860 


8(4) 


4670 (1200) 


5 (3) 


3860 (1200) 


2 (1) 


4180 (1000) 


3 (2) 


4700 (2000) 


3 (2) 


5700 (2500) 


Silll] 1892 


7(3) 


1900 (800) 


8 (2) 


1740 (700) 


6 (2) 


2530 (600) 


5 (3) 


2530 (1000) 


3 (2) 


3010 (1500) 


CIII] 1909 


10 (4) 


2350 (800) 


15 (4) 


1720 (700) 


17 (4) 


2510 (500) 


20 (7) 


3140 (700) 


18 (9) 


3920 (1200) 


Mgll 2798 


12 (8) 


2030 (700) 


14 (4) 


1500 (700) 


16 (2) 


1600 (600) 


40 (8) 


3960 (600) 


44 (8) 


6540 (600) 


Fcii uv 


22 (10) 


3000 (-) 


16 (8) 


3000 (-) 


9(5) 


3000 (-) 


7(5) 


5000 (-) 


7(5) 


9000 (-) 
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Table 3. CivA1549 components 



S.T. Component 1 Component 2 





Type 


A 


FWHM 


Strength (%) 


Type 


A 


FWHM 


Strength (%) 


A3 


Lor 


1547.3 


2220 


50 


Blue 


1532 


4800 


50 


A2 


Lor 


1547.8 


2600 


76 


Blue 


1534 


3900 


24 


Al 


Lor 


1547.7 


3580 


100 










Bl 


Gau 


1548.3 


4170 


42 


Gau 


1555.8 


15590 


58 


B1+ 


Gau 


1547.8 


4780 


28 


Gau 


1552.6 


17480 


72 
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Composite spectra for different bins 

1 2 I — I i n i 1 1 1 1 1 1 1 1 1 1 1 r~ 




Q I I I I [ U 1 I I I U I I U I I I I I I I I 1 I 

1000 1500 2000 2500 3000 

X Rest (Angstroms) 



Fig. 1. — Median HST spectra for the 5 spectral types considered in this paper, with iden- 
tification of the strongest emission lines. Each spectrum was corrected for the interstellar 
extinction and shifted to the rest frame before being processed. Horizontal scale is rest frame 
wavelength in Angstroms, vertical scale is in arbitrary F\ intensity units. Spectra have been 
offset on the vertical axis for clarity. 
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Composite spectra for different bins 



12 




1200 1400 1600 1800 2000 

X Rest (Angstroms) 

Fig. 2. — The rest frame region of the median spectra between 1150 A and 1950 A is shown 
after continuum and Felluv emission subtraction. As for Fig. 1, spectra have been offset for 
clarity. All lines for which we report equivalent width and FWHM in Table 2 are identified. 
Horizontal scale is rest frame wavelength in Angstroms. 
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1150 1200 1250 1850 1900 1950 

\est (Angstroms) 



Fig. 3. — Decomposition of blended lines for the spectral types considered in this paper. 
Left panels: Lja +NVA1240; right panels Al IIIA1860, Si III]A1892 and Cm]A1909. To all 
blends we subtracted in advance the underlying continuum and an Felluv emission template; 
Felluv emission is however weak around Lja and Cin]A1909. Thick lines show the profiles 
we fit; thin lines show individual line contributions (assumed to have Lorentzian profiles) to 
the blend. No narrow component is subtracted prior to the fitting in any of the lines. While 
the NVA1240 equivalent width remains approximately constant, the Si III] A1892/ Cin]A1909 
intensity ratio decreases monotonically along the sequence of spectral types. Horizontal scale 
is rest frame wavelength in A; vertical scale is intensity in arbitrary units. 
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CIV 1549 profiles for different bins 




1500 1550 1600 

X Rest (Angstroms) 



Fig. 4. — CivA1549 line profile decomposition for the spectral types considered in this 
paper. The fitted data are shown by thick lines; they are obtained from the continuum and 
Felluv subtracted spectra after He IIA1640 and Oin]A1663 and NIV]A1486 removal . For 
bin A3 and A2, the thin line superimposed to the data shows the best fit provided by an 
unshifted Lorentzian line to the CivA1549 red wing. The residual is a broad and blueshifted 
feature, possibly associated with outflows. For bin Al, a symmetric Lorentzian provides a 
satisfactory fit. For spectral types Bl and B1+, two Gaussian profiles, similar to the one 
used by Sulentic et al. in Paper I to fit H/^bc provide a satisfactory fits to the CIVA1549 
broad component. Narrow central residuals seen in Al and B-bins can possibly be associated 
with the contribution of the NLR. 
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RL - RQ difference for Lya and CIV (bin Bl) 

1 1 I 1 1 1 1 I 1 1 1 'I L 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 




1150 1200 1250 1500 1550 1600 1650 1700 

Kesi (Angstroms) 

Fig. 5. — Lya (left panels) and CivA1549 profiles (right) for radio-loud (top panels) and 
radio-quiet (middle panels) quasars. Note that, unlike in Fig. 3, contaminating lines have 
not been removed. Residuals (bottom panels) show that the main difference is due to a 
stronger narrow component in radio-loud sources. Horizontal scale is rest frame wavelength 
in Angstroms; vertical scale is in arbitrary intensity units. 
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Difference between brighter and faiter objects from bin Bl 




1200 1400 1600 1800 2000 

X Rest (Angstroms) 



Fig. 6. — Similar division as the one shown in Fig. 5 is presented here - this time we show the 
difference between the fainter objects and the brighter objects from bin Bl. The separation 
value of the absolute V-band luminosity My = —25 is chosen in order to create two equally 
represented groups. The Baldwin effect (anticorrelation of the equivalent width of a line with 
the luminosity is apparent for lines like Lya, Si IVA1398, CivA1549, Si III]A1892, but not 
for lines like NvA1240and Si II]A1306. The situation is unclear for He IIA1640, Olll]A1663 
and other lines due to their lower contrast. 



